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I develop a theoretical framework for inferring nonequilibrinm equations of motion from incom¬ 
plete experimental data. I focus on genuinely irreversible, Markovian processes, for which the 
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evolution, yet without any information about the timing of these snapshots. A reconstruction of 
the equation of motion must therefore be preceded by a reconstruction of time. 
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I. INTRODUCTION 

Markovian processes, both reversible and irreversible, 
are ubiquitous. Microscopic processes, described by the 
Schrodinger equation, are obviously Markovian; and so 
are many macroscopic processes that can be described 
by, say, a rate, master, or Boltzmann equation. All these 
processes share the common feature that they are local 
in time: The state of the system at any given time fully 
determines its future evolution, regardless of the system’s 
prior history; the system exhibits no memory. This ren¬ 
ders Markovian processes reproducible, in the sense that 
preparation of the same initial state — no matter how — 
always entails the same subsequent evolution. The ubiq¬ 
uity of Markovian processes is linked to the existence of 
disparate time scales in many systems. The macroscopic 
observables whose dynamics one wishes to describe typi¬ 
cally coincide with the slow observables, and hence they 
evolve on longer time scales; whereas any memory, which 
is due to interaction with the other, faster degrees of free¬ 
dom, fades away on a much shorter time scale. There are 
situations in which this separation of time scales breaks 
down and memory effects do play a role [1]; but in the 
present paper, I focus on Markovian dynamics. 

Much effort has been devoted to developing theoreti¬ 
cal frameworks that allow one to derive Markovian trans¬ 
port equations from the underlying microscopic dynam¬ 
ics; and in turn, to deduce from these transport equations 
testable predictions for macroscopic experiments Hi]. 
In contrast to this deductive approach, I start here from 
the opposite end: I ask how one can infer a Markovian 
transport equation/rom experimental data. In particular, 
I consider situations where the experimental data come 
without time information. Such data may stem from past 
processes (say, in the geological or astronomical realm) 
that have left visible traces, albeit without time informa¬ 
tion; or from processes that, again, leave visible traces 
but are so fast or delicate that they cannot be tracked 
with a clock. Under such circumstances, before inferring 
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the pertinent Markovian equation of motion, one must 
first reconstruct “time.” The purpose of the present pa¬ 
per is to show that this is possible (up to an additive 
and multiplicative constant); and to furnish the neces¬ 
sary tools for doing so. 

The reconstruction of time and of the equation of mo¬ 
tion presupposes rather detailed knowledge about the 
generic structure of Markovian dynamics. In particular, 
Markovian dynamics may be viewed geometrically as a 
flow on the manifold of macrostates. This manifold is 
endowed with a rich geometric structure, and there are 
certain consistency conditions that any Markovian flow 
must satisfy. These conditions, in conjunction with the 
experimental data, will turn out to constrain the form 
of the equation of motion just enough so that the recon¬ 
struction succeeds. Mirroring the importance of these a 
priori constraints, I start out with a comprehensive dis¬ 
cussion of the geometry of macros tate s (Sec. 0 and of 
generic Markovian dynamics (Sec. HI)- These introduc¬ 
tory sections draw on ideas from the projection-operator 
HD, geometric [7 ]j two-generator mm approaches 
to nonequilibrium dynamics. Then I proceed to formu¬ 
late the general prescription for reconstruction (Sec. IVI, 
which I later illustrate with a simple example (Sec.^l. 
I conclude with a brief discussion and outlook on future 


work. Sec. VI 


II. GEOMETRY OF MACROSTATES 
A. Manifold and coordinates 

To describe the static or dynamic properties of a 
macroscopic system, typically only a few observables 
are deemed relevant - for example, the system’s con¬ 
stants of the motion (if static), slow observables (if dy¬ 
namic), or observables pertaining to some subsystem of 
interest. These relevant observables, together with the 
unit operator, span the so-called level of description, 
a subspace within the linear space of observables [Hj. 
For an arbitrary microstate p and level of description 
Q := span{/, Ga}, the associated macrostate is that state 
which, while yielding the same expectation values {ga} 
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as p for all relevant observables, comes closest to equidis- 
tribution. Closeness to equidistribution is measured in 
terms of the von Neumann entropy 

S[p] :=-tiiplnp); (1) 

so the macrostate - denoted by 7r(p) - is determined by 
the maximization 


7r(p) := argmaxS'[/x], (2) 

ueg 

where p € g is short for the constraints {Ga)^ = goSo,. It 
has the Gibbs form 

7r(p) = ^(A)“^ exp(-A“Go), (3) 

with the partition function 

Z(X) := tr{exp(-A“Ga)} (4) 


which in turn implies that the Hessian of the entropy 
yields (up to a sign) the inverse of the correlation matrix, 

((^-l)ab ^ _d‘^dbg (g) 

Expectation values x := (Af),n.(p) of arbitrary (not neces¬ 
sarily relevant) observables X change by 

dx = dx{d^)dgb = -{SGb;X)dX\ (10) 


where 6Gb '■= Gb — gb, and (;) denotes the canonical 
correlation function 


{A;B) := f dj^tr[7r(p)^HV(p)^ B], (11) 

Jo 

The latter constitutes a positive definite scalar product in 


the space of observables. Comparing Eq. (10) for a; = 
with Eq. (|^, one finds that the correlation matrix can 
be expressed in terms of this scalar product. 


ensuring state normalisation, the Lagrange parameters 
{A“} adjusted such that {Ga)-K(p) = ffa, and - for ease of 
notation - the Einstein convention that identical upper 
and lower indices are to be summed over. The opera¬ 
tional meaning of “relevance” and how it leads to the 
Gibbs form have been discussed in Ref. m- 

Let S denote the set of normalized (pure or mixed) 
states of a given physical system; this set constitutes a 
differentiable manifold. In this manifold, the macrostates 
with level of description Q form a submanifold, 7r(iS). 
This submanifold has dimension (dim^—1), which equals 
the number of relevant observables (provided they are 
linearly independent). On 7r(5) there are two natural 
choices of coordinates: the relevant expectation values 
{9a ■■= {Ga) ii-(p)}, or the Lagrange parameters {A“}. The 
former coordinates can be expressed in terms of the latter 
via 


Gab = {5Ga-,5Gb). ( 12 ) 

B. Coarse graining and projectors 

The map tt : 5 —^ ^(>5) constitutes a coarse graining. 
While retaining information about the relevant observ¬ 
ables, it discards all information about the rest. Geo¬ 
metrically, it “projects” a microstate p G S onto the sub¬ 
manifold of macrostates. Indeed, tt exhibits typical fea¬ 
tures of a projection operator: it is idempotent, ttott = tt; 
successive coarse grainings with respect to smaller and 
smaller levels of description are equivalent to a one-step 
coarse graining with respect to the smallest level of de¬ 
scription, 

G CB ^ TTg OTTjr = TTg] (13) 


ga = -da In Z, (5) 

where da '■= d/dX°'. Their respective gradients {dga} and 
{dX°’} are related by 

dga = -GabdX^ , dA“ = -{G-^r^dgb, (6) 

with the Jacobian (up to a sign) given by the correlation 
matrix 

Gab-dadblnZ. (7) 

Associated with the expectation value coordinates {ga} 
is a local basis of tangent vectors := d/dga], related 
to their one-form duals {dga} via dgb{d°‘) = 5^ [TT]; and 
likewise for the Lagrange parameter coordinates, with the 
associated local basis {9o} satisfying dX^{da) = 6{^. 

Upon infinitesimal variation of a macrostate its von 
Neumann entropy changes by 

( 8 ) 


and it is covariant under unitary transformations, 

7Tuguf{UpU^) = UTTg{p)UK (14) 

In contrast to an ordinary projection operator, however, 
this coarse graining map need not be linear. 

Mirroring the coarse graining of states, there is also 
a coarse graining (“super-” )operator V on the space of 
observables. The two are dual to each other in the sense 
that for arbitrary (not necessarily relevant) observables 
X and arbitrary microstates p it is 

{vx)p = {x)^^p) yx,p. (15) 

The unique superoperator which satisfies this require¬ 
ment, sometimes called the Robertson |5] or Kawasaki- 
Gunton [12j projector, is 

rx ■.= xI + dx{d'^)6Ga. (16) 

It projects arbitrary observables onto the level of de¬ 
scription, the projection being orthogonal with respect 


dS = dS{d^)dga = X^dga; 
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system 

discard 

7r(p) 

vx 

single 

coherence 


EiP^XB 

SxE 

environment 

PS ® bE 

tlE(bEX) g) Ie 

Ax B 

correlations 

pA g) PB 

tXBiPsX) g) Ib 
+Ia g) tr a{paX) 
-xIa g) Ib 


TABLE I. Three examples of levels of description and the 
associated coarse graining operations tt and V. They refer 
to (i) a single quantum system where only classical proba¬ 
bilities are deemed relevant, and all information about co¬ 
herence is discarded; (ii) a system S coupled to an environ¬ 
ment E where only the system properties are deemed rel¬ 
evant, and all information about the environment (includ¬ 
ing system-environment correlations) is discarded; and (iii) 
a bipartite system where only single-particle properties are 
deemed relevant, and all correlations are discarded. The {P;} 
are projection operators (on Hilbert space) pertaining to some 
preferred orthonormal (“decoherence”) basis, pi denotes the 
reduced state of subsystem i. The state t := I/txI is the to¬ 
tally mixed state. In the hrst two examples tt is linear, and 
V is state-independent. In contrast, in the third example tt is 
nonlinear, and V varies with the macrostate. 


to the scalar product (;). As the latter is evaluated in 
the macrostate 7r(p), and hence the notion of orthogonal¬ 
ity varies with the macrostate, the projector, too, carries 
an implicit dependence on the macrostate: V = 'P[7r(p)]. 
Like the coarse graining operation tt on states, the projec¬ 
tor is idempotent, = V; yet unlike tt, it is always linear 
and hence a true projector. Its complement Q := 2 — V 
(with I being the unit superoperator, IX = A) is also 
a projector and projects an arbitrary observable onto its 
“irrelevant” component. Table |T] lists some examples for 
different levels of description and the associated coarse 
graining operations tt and V. 


C. Metric and covariant derivative 


As the canonical correlation function (;) is a positive 
definite scalar product, Eq. (121 implies that the cor¬ 
relation matrix is both symmetric and positive definite. 
Therefore, the symmetric (0,2) tensor field 

C := Cab (g) dX^ = (C-i)“'’d5a g) dgt (17) 


constitutes a Riemannian metric on the manifold of 
macrostates. This particular metric is known as the Bo- 
goliubov or Kubo-Mori metric [laiii]. Up to a sign, it 
relates the local basis in ^-coordinates to the one-form 
duals in A-coordinates and vice versa, 

C(a“) = -dA“ , C{da) = -dga. (18) 

The Bogoliubov-Kubo-Mori (BKM) metric has an op¬ 
erational meaning. It quantifies the statistical distin- 
guishability of nearby macrostates, in the following sense. 


Two states p and p can be distinguished statistically if 
measurements on a finite sample, taken from an i.i.d. 
source of one state, say, p, are highly unlikely to erro¬ 
neously indicate the other state, p. The pertinent error 
probability is 

probi_,(^|A, p) 

:= inf {prob(r|p®^) | prob(r|^®^) > 1 - e} , (19) 

where N denotes the size of the sample, and T is a proxy 
for the measurement results, which asymptotically, i.e., 
to within an error probability e (0 < e < 1) that does 
not depend on sample size, are compatible with the sam¬ 
ple being in the state p'^^ . Asymptotically, this error 
probability decreases exponentially with sample size, 

probi_,(/r|iV,p) - exp[-AS'(/r||p)], (20) 

and no longer depends on the specific value of the er¬ 
ror parameter e (“quantum Stein lemma” [13 HD). The 
exponent features the relative entropy iniiii], 

II , I tr(/iln/i — /rlnp) : supp/i C supp p 

^ ^ [ -boo : otherwise 

( 21 ) 

which thus proves to be a natural distinguishability mea¬ 
sure. For nearby macrostates 7r(p) and 7r(p -|- 6p) con¬ 
nected by a distance vector W, this distinguishability 
measure is approximated to lowest (quadratic) order in 
the coordinate differentials by 

S{nip)Mp + Sp))^lC{W,Wy, (22) 

and hence indeed, up to a numerical factor, by the length 
(squared) of the distance vector in the BKM metric. 

The above metric is also singled out by the fact that 
it is with respect to this metric that the “projection” tt 
onto the submanifold of Gibbs states is orthogonal. This 
can be seen as follows. Being the infinitesimal version 
of the relative entropy, and the latter being defined for 
arbitrary pairs of states (both macro and micro), the 
metric can be extended from the manifold of macrostates 
to full state space. Let V denote an arbitrary vector field 
on the full state space that connects only states with 
identical expectation values for the relevant observables, 
dga{y) = 0. Then with dga = —C{da), it is 

C(aa,U)=0Va; (23) 

i.e., any such V intersects 7r(5), which is generated by 
the basis vectors {da}, at a right angle. 

A Riemannian metric allows one to raise and lower in¬ 
dices of tensor fields, i.e., to map an (n, m) tensor field to 
an (n-|-l,TO—1) or (n—l,m-|-l) tensor field, respectively, 
thereby preserving its total rank, n+m. Henceforth I will 
not distinguish between tensor fields that differ only by 
raising or lowering of indices via the BKM metric; I will 
denote such fields by the same symbol and characterize 
them only by their total rank. 
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Associated with the metric is a covariant derivative, 
V. Given some vector field V, the covariant derivative 
along V, Vy, maps an arbitrary rank-r tensor field A to 
another rank-r tensor field, VyA. This map (i) obeys the 
sum rule, Vy(A + B) — (VyA) -f (VyS); (ii) obeys the 
Leibniz rule for tensor products, Vy(A(g) B) = (VyA) (g) 
B + (VyB); (iii) commutes with tensor contraction; 

(iv) is linear in the vector field, V/[/+gyA = fWuA + 
gVyA, for arbitrary scalar functions /, g and vector helds 
U, V. For this reason there exists a rank-(r -|- 1) tensor 
field, denoted VA and called the gradient of A, from 
which VyA can be obtained via contraction with V; and 

(v) when applied to a scalar field 4>, the map coincides 
with the ordinary directional derivative, Vy = dct>{V). 
So in this special case the gradient coincides with the 
ordinary differential, = dcj). 

That the covariant derivative stems from the BKM 
metric is reflected in the fact that the gradient of the 
metric tensor vanishes, 


VC = 0. (24) 

In particular, there is no torsion. 


VuV - VyU = [U, V] V U, V, (25) 

the bracket [U, V] being the Lie bracket of the two vector 
fields U and V [U. The absence of torsion implies that 
the gradient of a one-form field a can be written as 

Va= ^[da + £c-iic.)C], (26) 

where d denotes the exterior derivative and £ the Lie 
derivative. The first term inside the square bracket is 
antisymmetric, whereas the second term is symmetric. 
Whenever a is itself a gradient, a = dcj) = V(j), its exterior 
derivative vanishes, and hence it is 

VVcl)=\£c-ny4,)C. (27) 

One scalar function that will play a special role in my 
subsequent argument is the modihed entropy [19] 

S„[g] := S[a]-S{^i\\a) (28) 

with reference macrostate a € This modified 

entropy characterizes the closeness of /r to the refer¬ 
ence macrostate; it reduces to the ordinary von Neu¬ 
mann entropy when the reference macrostate equals the 
totally mixed state. Upon infinitesimal variation of a 
macrostate (at fixed reference macrostate) the modified 
entropy changes in a manner similar to the ordinary en¬ 
tropy, Eq. 0, only with A“ replaced by (A“ — A“), 

dB, = (A“ - A“)d5„ (29) 


where the {A“} pertain to the reference macrostate a. 
Taking the gradient of this differential yields, with the 
help of Eqs. (§ and 0, 


VVB, = - 


(30) 


So when evaluated at the reference macrostate, where 
A“ = A“, this covariant Hessian equals (up to a sign) the 
metric tensor. 


VVS,y = -C at cr. (31) 


When macrostates are constrained to the hyperplane 
T, := {fi G 7r(B) I (In cr)^ = (Incr)cr}, modified and ordi¬ 
nary entropies coincide, Bo-js = B|e. Then Eq. (31) 
carries over to the ordinary entropy. 


VVB|e = -GIe at cr. 


(32) 


As an example, a might be a canonical equilibrium state 
and hence In cr proportional to the Hamiltonian. Then S 
constitutes a hyperplane of constant energy. It comprises 
all macrostates that have the same energy as a, includ¬ 
ing cr itself. Whenever energy is conserved, evolution of 
the macrostate is constrained to such a hyperplane; so 
energy-conserving flows have C\y, as their relevant met¬ 
ric. At the equilibrium state cr this constrained metric is 
given by the (negative) Hessian of ordinary entropy. 


III. GENERIC MARKOVIAN DYNAMICS 
A. Disparate time scales 

In this paper I focus on the dynamics of an isolated 
quantum system with time-independent Hamiltonian H. 
(In principle, the description of an open system can be in¬ 
corporated into this framework by enlarging it to include 
its environment.) On the microscopic level the dynamics 
of such a system is governed by the Liouville-von Neu¬ 
mann equation 


p{t) = -iCp{t), (33) 

where p denotes the system’s microstate, and the Liouvil- 
lian C := [H, •] is a shorthand for the commutator with 
H. For simplicity, I set h = 1. 

On the macroscopic level one seeks to describe the 
dynamics of only certain selected expectation values 
ga{t) := (Ga)p(t) deemed “relevant”. Provided that ini¬ 
tially, at t = 0, these relevant expectation values suffice 
to determine the system’s microstate, i.e., p(0) carries no 
information other than about {5a(0)} and hence has the 
Gibbs form 


p{0) cx exp(-A“(0)Ga), (34) 

their dynamics at t > 0 is governed by the Robertson 
equation 


9a{t)=gi^\t)+gi^\t) (35) 

with the local term 

9a\t) = 


G + (A“ - A“)VVg, 


(36) 
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and the memory term 

g^rHt) = - f dt'{CQ{t')nt\t)Q{t)CGa).(pit')). 

(37) 

Here Tr{p{t)) is the macrostate at time t as defined in Sec. 
|II A[ The objects Q and T are, like the Liouvillian, super¬ 
operators acting on the space of observables. The former 
is a projector, = Q, which projects any observable 
onto its “irrelevant” component; it is the complement of 
the Robertson projector defined in Sec. m Like the 
Robertson projector, it may vary with the macrostate 
and thus may carry an implicit time dependence. The 
superoperator T effects the time evolution of irrelevant 
degrees of freedom, 

id/dt')r{t', t) = -iQ{t')CQ(t')T{t', t), (38) 


with initial condition Tit, t) = I. 

All terms on the right-hand side of the Robertson equa¬ 
tion depend on the macrostate and hence on relevant ex¬ 
pectation values only; so the system of equations of mo¬ 
tion for the {5'a(i)} is indeed closed. Irrelevant degrees 
of freedom have been eliminated completely from the de¬ 
scription of the macroscopic dynamics. The price to pay 
for this elimination is that in contrast to the microscopic 
Liouville-von Neumann equation, the Robertson equa¬ 
tion can be both nonlocal in time and nonlinear. The for¬ 
mer means that the change of relevant expectation values 
at any given time may depend not just on their current 
values but on their entire history since t = 0, i.e., that 
the macroscopic dynamics has a “memory”. The latter 
- nonlinearity - may arise whenever the coarse graining 
operation tt is not linear. A simple example for such a 
nonlinear coarse graining was given in Table |l] Indeed, 
many well-known transport equations such as the Boltz¬ 
mann or Navier-Stokes equations, which can be derived 
within the above framework or its classical counterpart, 
are nonlinear. 

The Robertson equation becomes Markovian, i.e., local 
in time, if and only if the physical system exhibits a clear 
separation of time scales, and it is the slow degrees of 
freedom which are chosen as the relevant ones. In this 
case the “memory time” - the time scale on which the 
integrand in Eq. (371 falls off towards the past - will be 
short compared to the time scale r^. on which the relevant 
expectation values evolve. One may then replace 


where for simplicity I omitted the dependence on t. 

Geometrically, the collected relevant expectation val¬ 
ues at any given time can be represented as a point, 
and their time evolution as a curve, in the manifold 
of macrostates. The curve results from projecting the 
trajectory of the microstate in full state space onto the 
lower-dimensional submanifold of macrostates. In case 
the Robertson equation is Markovian (which I shall as¬ 
sume from now on), its general solution defines in the 
manifold of macrostates a congruence of curves. This 
congruence in turn gives rise to a vector field on the man¬ 
ifold of macrostates. 


V-.= gad\ 


(41) 


The correspondence being one-to-one, the Markovian dy¬ 
namics may be characterized completely by the vector 
field V. In line with the split in Eq. (35), the vector field 
can be broken down into contributions from the local 
term and the memory term, V = 


B. Effective non-dissipative dynamics 

The effective non-dissipative dynamics of the relevant 
expectation values is described by the local term and the 
antisymmetric part of the memory term. As for the local 
term, one can use the general formula 

[X,p\= f di,p''[X,\np\p^-’' (42) 

Jo 

and define on the manifold of macrostates the antisym¬ 
metric (2,0) tensor field 

K := C-[Ga, G,]; {H - (i7).(p))).(p)5' G 9“ (43) 

to cast it into the form 

=K{dS,-). (44) 

As an immediate consequence of the antisymmetry of K, 
the local term preserves entropy, 

K{dS, dS) = 0, (45) 

and is thus indeed non-dissipative. 

As for the memory term, using the relation 


T^{p{t)) > Q[i') Q(t) ( 39 ) 

in both the Robertson equation and the differential equa¬ 
tion for T (“Markovian approximation”). Moreover, pro¬ 
vided there is a genuine gap between the two time scales, 
in the sense that there exists an intermediate scale T such 
that Tjn Tr, and with the substitution {t — t') —>■ t 

one may expand the integration range for r from [0, t] to 
[0, T] and replace 

f dt'T{t',t) ■= f dr exp{iTQCQ), (40) 

Jo Jo 


{CX)^(p)=-\\CGb-X)^(p) (46) 

for arbitrary (not necessarily relevant) X, the hermiticity 
of Q with respect to the canonical correlation function, 
as well as the fact that the memory term must be real, 
and defining the (2, 0) tensor field 

M := (l(+) QCGa, QCGb)^(^p)d^ ® (47) 

it can be cast into the form 

rW = M{dS, ■). 


(48) 
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The tensor M comprises an antisymmetric and a sym¬ 
metric part whose respective components are given by 

(49) 

The antisymmetric part conserves entropy, 

M^-'>{dS,dS) =0, (50) 

and thus co-determines, together with the local term, the 
non-dissipative dynamics. Often neglected, this antisym¬ 
metric part of the memory term may contain interesting 
physics such as geometric or other effective forces m, 
and it may play an important role in Hamiltonian renor¬ 
malization |21| . The most general Markovian dynamics 
is then described by the tensor 

T := {K + M^-^) + , (51) 

where the first two, antisymmetric terms drive the non- 
dissipative dynamics, whereas the last, symmetric term 
drives the dissipative dynamics. With this tensor the 
complete equation of motion acquires the compact form 

V = T(dS, •)• (52) 


It is possible to describe the non-dissipative dynamics 
by an effective Hamiltonian if and only if 


(K + M(-^)(dS,-)=K,ff(dS,-), 


(53) 


where Keg is defined as in Eq. ( |43[ ) but with H replaced 
by Heg. In this case the non-dissipative time evolution 
of macrostates is unitary. Since unitary transformations 
leave relative entropies invariant and hence, thanks to 
Eq. (22), also the BKM metric, it is 


£ 


K^tddS,-) 


.)C = 0. 


(54) 


Provided Keff is itself a relevant observable, the non- 
dissipative dynamics can be written in the alternative 
form 


Keg(dS,-)=L(d[/,-), (55) 

where U := (i?eff) 7 r(p) is the (effective) internal energy, 
and L denotes another antisymmetric (2, 0) tensor field, 

L:= (jlGa,Gh]U(,)d'^^d£ (56) 

In this formulation the conservation of entropy is re¬ 
flected in the tensor property 


L(-,dS) = 0. (57) 

In principle, the effective Hamilton operator may be in¬ 
fluenced by the macrostate, ffeff = and thus, 

through the latter, depend on time. If, however, such an 
explicit time dependence is absent or may be neglected, 
the antisymmetry of L ensures the conservation of inter¬ 
nal energy. 


(58) 


In case the relevant observables form a Lie algebra, 
(l/i)[^,0] C G, the manifold of macrostates and their 
non-dissipative dynamics exhibit further structure. Such 
a Lie algebra property holds for many important choices 
of the level of description: for instance, when the relevant 
observables comprise (i) all constants of the motion; (ii) 
all observables pertaining to one or several subsystems 
of a composite system; or (iii) all block diagonal observ¬ 
ables of the form '^■PiAPi, where {Pi} is some set of 
mutually orthogonal projectors. In all these examples 
the commutator (times l/i) of two relevant observables 
is again a relevant observable. Then the bracket 

{f, 9 } ■= L{df,dg) (59) 

of two functions / and g, defined with the help of the ten¬ 
sor field L, possesses all properties of a Poisson bracket; 
it satisfies antisymmetry, linearity, Leibniz rule, and the 
Jacobi identity. Thus the manifold of macrostates is en¬ 
dowed with a Poisson structure [22] . According to the 
splitting theorem for Poisson manifolds [23] the mani¬ 
fold of macrostates can then be foliated into symplectic 
leaves^ each with constant entropy. On every leaf one can 
define a symplectic two-form, i.e., a two-form which is an¬ 
tisymmetric, non-degenerate, and closed. Also on every 
leaf, the vector field associated with the non-dissipative 
dynamics, L(dU,-), becomes a Hamiltonian vector field, 
the pertinent Hamilton function being the internal energy 
U. One recovers thus the familiar structure of classical 
Hamiltonian mechanics [24] . 


C. Dissipation 


Dissipation is described by the symmetric part of the 
memory term. This symmetric part is positive semidefi- 
nite, 

> 0, (60) 


which can be understood as follows. One assumes that 
on short time scales smaller than T the dynamics of the 
irrelevant degrees of freedom is (i) time translation in¬ 
variant; and in particular, (ii) unaffected by the (slow) 
variation of the macrostate. Then it is 


(Q£G,; J(+)Q£G6) = [I^-^QCGa] QCGb), (61) 


where is defined as in Eq. (40) but with a minus sign 
in the exponent. Considering moreover that the memory 
term must be real, and hence that the canonical correla¬ 
tion function (which is a scalar product) featuring in the 
definition of M must be symmetric, the components of 
the symmetric and antisymmetric parts of M are given 

by 

^£^-^)QCGa;QCGb). (62) 


Eor the symmetric part one can then write (invoking once 
more the time translation invariance of the irrelevant dy- 


L{dU, dU) = 0. 
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namics on short time scales) 


M[t'> = 


dr {exp{tTQCQ)QCGa; Q£Gb) 


f-T 




dr 


-T 


x(exp[i(T + s)QCQ]QCGa] ex-p{isQLQ)QLGb) 


pT+s 


ds 


dr' 


l-T+s 


2T . 

X (exp (ir' Q£ Q) Q/:G, ; exp(zs QC Q) QCGb). 


(63) 


To the second integral only t' s G [0,T] contribute 
significantly, and hence one may reduce its integration 
range [—T + s,T + s] —)• [0, T]. This finally yields 


= (l/2r)(l(+)Q£G,;l(+)Q/:Gb). (64) 


Since the canonical correlation function is a positive def¬ 
inite scalar product, this component matrix, and hence 
itself, must indeed be positive semidefinite. 

An immediate consequence of the above positivity is 
that the symmetric part of the memory term may only 
lead to an increase, but never to a decrease, of entropy, 

M^+'>{dS,dS)>0. (65) 


This embodies the H-theorem, originally formulated by 
Boltzmann for the dynamics of dilute gases [21] but in 
fact valid for arbitrary Markovian processes. It marks the 
gradual approach of the macrostate towards equilibrium. 

One particularly simple type of dissipative dynamics 
is steepest descent towards the equilibrium state. Under 
this dynamics trajectories on the manifold of macrostates 
(or in case there are conservation laws: on the allowed 
submanifold) simply follow the entropy gradient |26j . 
The only contribution to this type of dynamics stems 
from the symmetric part of the memory term, which has 
the simple form 

oc G-\ (66) 


This is a rather natural ansatz: The metric G measures 
the distance between two macrostates in the sense of their 
statistical distinguishability. The symmetric part of the 
memory term (or its inverse, respectively), being sym¬ 
metric and positive semidefinite, is a metric, too; it mea¬ 
sures the “dynamical” distance between two macrostates 
as mediated by interactions with the irrelevant degrees 
of freedom. In the absence of any specific information 
about the dynamics of the irrelevant degrees of freedom, 
in particular about any preferred direction on the man¬ 
ifold of macrostates, these two metrics are taken to be 
equal, up to some multiplicative constant. A direct con¬ 
sequence of this ansatz is that the gradient of van¬ 
ishes everywhere, = 0. This is not the case for 

other, more general forms of Markovian dynamics. The 
magnitude of this gradient can then be taken as a local. 


coordinate-independent measure for the deviation from 
steepest descent. 

When the overall dynamics drives the macrostate to¬ 
wards an equilibrium macrostate a € Tr(S) then its log¬ 
arithm, In cr, which is a relevant observable, must be a 
constant of the motion. Hence in their approach towards 
CT, macrostates are constrained to the hyperplane E de¬ 


fined in Sec. IIC Provided the equilibrium macrostate 
is a canonical state with the same effective Hamilto¬ 
nian as that governing the non-dissipative dynamics, 
a oc exp(—/3a-Jdes), this hyperplane corresponds to fixed 
internal energy, [/ = const. The latter is conserved by 
the non-dissipative dynamics, Eq. (58), and so in order 


to be conserved overall, it must be conserved by the dis¬ 
sipative dynamics, too. This conservation of energy is 
implemented mathematically by imposing 


M^+\-,dU) = 0. 


(67) 


D. Near equilibrium 


Let cr G TT (5) be an equilibrium macrostate a nd S C 

to 


7r(5) the associated hyperplane, defined in Sec. HC 


which macrostates are constrained as they approach a. 
In the following, all mathematical objects —macrostates, 
functions, vector and tensor fields— and their relation¬ 
ships are meant to be constrained to this hyperplane, 
even if for simplicity I will omit the explicit notation 
“Id”. To begin with, taking the covariant derivative on 
both sides of the general equation of motion, Eq. (521, 


and exploiting the relationship between the Hessian of 
the entropy and the metric tensor, Eq. (32), as well as 


dS = 0 at equilibrium, yields an equation for T at equi¬ 
librium. 


T = -VU at ct; (68) 

i.e., at equilibrium, the tensor T which governs the 
Markovian dynamics can be gleaned from the observed 
vector field V by taking (minus) the gradient. To first- 
order approximation, this relationship can be extended 
to the vicinity of equilibrium, 

T « -VU; (69) 

which in turn, by the general equation of motion, implies 

U«-(VU)(dA,-). (70) 


The latter equation imposes a consistency condition on 
the vector field V near equilibrium. It goes beyond the 
trivial requirement that U = 0 at u. 

By symmetrizing both sides of Eq. (69) and exploiting 
the positivity of the symmetric part of T, Eq. (601, one 
finds that the symmetric part of the gradient VU must 
be negative semidefinite. This symmetric part is denoted 
by [VU]+ and may be called, in an analogy with fluid 









dynamics, the rate of strain tensor. In conjunction with 
Eq. (26), its negativity leads to the inequality 


£vC = 2 [VE]+ < 0, 


(71) 


again valid in the vicinity of equilibrium. Near equilib¬ 
rium, therefore, macrostates converge (in terms of the 
BKM metric) not just towards equilibrium but also to¬ 
wards each other; as they evolve, their mutual distances, 
and hence their statistical distinguishability, decrease 
monotonically. This inequality may be viewed as a spe¬ 
cial case of Lindblad’s theorem m, here applied to the 
linearized dynamics near equilibrium. 

When the dynamics is entirely non-dissipative, = 
0, the above reasoning implies £vC = 0. This is in 


keeping with the unitarity condition, Eq. (541. As the 


BKM metric is the infinitesimal version of relative en¬ 
tropy, Eq. ( |^ , its conservation suggests that, more¬ 
over, arbitrary relative entropies are preserved. If this 
(unproven) conjecture were true then the conservation 
of the BKM metric, £vC = 0, would in fact mandate, 
rather than merely allow, that the non-dissipative dy¬ 
namics is unitary [351ISH] ■ 


IV. RECONSTRUCTION 

The general framework laid out above offers a sys¬ 
tematic route from the microscopic to the macroscopic 
realm. Starting from the microscopic Hamiltonian and 
the Liouville-von Neumann equation, Eq. (33), it al¬ 


lows one, at least in principle, to calculate the tensor 
field T that governs the Markovian macroscopic dynam¬ 
ics; and from there, via the general equation of motion, 
the vector field V, and hence the time evolution of macro¬ 
scopic observables. Reconstruction starts from the other 
end: Having observed some macroscopic Markovian pro¬ 
cess experimentally, one aims to infer the dynamical law. 
Specifically, one endeavours to reconstruct the pertinent 
tensor field T. 

The generic experimental setting can be characterized 
as follows. On some coarse-grained level of description, 
one observes a hitherto unknown process whose only 
known feature is that it is reproducible and hence Marko¬ 
vian. The dynamical law governing this macroscopic pro¬ 
cess - here: the tensor field T - is not known and yet 
to be inferred from the data. It is impossible to infer 
the dynamics from a single trajectory alone; rather, one 
has to study an entire family of trajectories pertaining 
to differently prepared test systems. More specifically, 
in order to collect sufficient data for such an inference, 
one must (i) observe multiple copies of the system of in¬ 
terest, or provided the system is large enough so that 
disturbances due to measurement may be neglected, ob¬ 
serve one system undergo the same process several times; 
(ii) vary the initial macrostates of the different copies or 
sequential runs, respectively; (iii) in each run, measure 
the macrostate at different (ideally, closely spaced) in¬ 
stances of its evolution; and (iv) record the time between 


subsequent measurements. In combination, these data 
allow one to determine the vector field U, as illustrated 
schematically in Fig. [^a). 

Having thus determined experimentally the vector field 
U, one proceeds to infer the macroscopic dynamical law, 
embodied in the tensor field T. Near equilibrium, this is 


a straightforward exercise: According to Eq. (69), T is 


but the negative gradient of V. Calculating the gradient 
requires knowledge of the (generally non-Euclidean) ge¬ 
ometry of the manifold of macrostates, which depends on 
the level of description but is otherwise independent of 
the experimental data. Further away from equilbrium, 
there may be corrections to Eq. (69). One will then 


end up with some system of first-order partial differen¬ 
tial equations for the components of T. The difficulty of 
solving these differential equations will vary depending 
on the system at hand. 

For simplicity, I stick here to the vicinity of equilibrium 
and focus instead on another issue that may complicate 
the reconstruction of T. I consider the situation where 
one is given data points on the manifold of macrostates 
but without any information about time (Fig. 0 b)). On 
a large scale, this might be the case for, say, geological 
or astronomical data originating from different epochs 
of terrestrial or cosmic evolution, respectively, that are 
as yet undated but presumed to be connected via some 
hitherto unknown Markovian process (say, some yet-to- 
be-discovered form of stellar evolution). Or at the oppo¬ 
site extreme, it might apply to processes on a very small 
scale (say, in particle or molecular physics) that leave a 
visible trace but are too fast or delicate to be tracked 
with a clock. More exotically, such timeless data might 
be collected in a remote lab without access to an exter¬ 
nal reference clock and with all internal clocks broken, 
so that the very notion of time has to be reconstructed 
from scratch. In all these cases it is clearly impossible 
to construct the vector field V, and hence to infer the 
tensor T, in the straightforward manner outlined above. 

Nevertheless, with the help of assumptions to be dis¬ 
cussed below, a reconstruction of the Markovian dynam¬ 
ics is feasible even under such adverse circumstances. 
One important prerequisite is that the dynamics is gen¬ 
uinely irreversible, 5" > 0. Then at least the temporal 
order of the snapshots is easily established; and the en¬ 
tropy being a strictly monotonic function of time, it may 
in a first step serve as a proxy for time. With the help of 
this proxy one can construct a vector field W, related to 
the true V hy V = SW (Fig. 0c)). The reconstruction 
task thus reduces to the problem of finding the local en¬ 
tropy production rate, S. In the following, this unknown 
entropy production rate shall be denoted by a separate 
letter, 77 . 

The true, yet to be determined vector field V must 
satisfy a consistency condition, Eq. ( |70| ). Replacing V 
by rjW in this condition and exploiting dS{W) = 1 leads 
to a first constraint on 77 , 


dr] ■ dS ~ —277, 


(72) 
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FIG. 1. Schematic illustration of the basic reconstruction 
idea. The figures show a section of the hyperplane E of 
macrostates, (a) Experimental data are given in the form 
of timed snapshots of the macrostate at different instances 
of the Markovian evolution. Crosses and dots represent data 
from two different runs of the experiment with varying initial 
conditions; in practice, there will be many more such runs 
and data points. In combination, these data allow one to 
construct the vector field V. (b) When time information is 
removed, it is no longer possible to construct the vector field 
V. (c) However, as long as the dynamics is genuinely irre¬ 
versible, 5 > 0, at least the temporal order of the snapshots 
is easily established, and one may use entropy as a proxy for 
time. Here concentric circles around the equilibrium state a 
indicate hypersurfaces of constant entropy. With this proxy 
one can construct a vector field W, related to the true V by 
V = SW. The entropy production rate, S, can in turn be 
estimated from the data, up to a multiplicative constant, if 
one assumes the existence of an effective Hamiltonian for the 
non-dissipative part of the dynamics (see text). 


where the scalar product dr] ■ dS is short for C~^{dr], dS). 
I presume that the non-dissipative part of the dynamics 
is governed by some effective Hamiltonian, and that at 
least in the vicinity of equilibrium, this effective Hamil¬ 
tonian does not vary with the macrostate. The macro¬ 
scopic non-dissipative dynamics thus shares with the mi¬ 
croscopic dynamics a common structure, namely, unitary 
evolution with a time-independent Hamiltonian. Such 
structural invariance is not guaranteed a priori; rather, 
it constitutes an extra assumption which, however, is jus¬ 
tified for many systems. Mathematically, by Eqs. ( [5^ 
and (69), and taking into account the change of sign when 
going from 9“ to dA“ via the metric, Eq. (18 1 , this as¬ 
sumption means that there must exist an iJeff such that 


[VE]_(d5, •) « (z£effG„),(p)dA“. (73) 

Here [VE]_ denotes the antisymmetrized gradient 
(“curl”) of V. Again replacing V by rjW, exploiting 
dS{W) = 1, and using the previous constraint on 77 leads 
then to the condition 


^dr] K,r]uj- {iCeffGa)Tr(p)dX°-. (74) 

On the right-hand side, the one-form 

oj ■.= [VW]-{dSr)-W (75) 

is determined by the observed proxy field W. The ef¬ 
fective Hamiltonian (up to an additive and multiplica¬ 
tive constant), and hence the effective Liouvillian and 
the expectation values {iLesGa) (up to a multiplicative 
constant), can be inferred from the observed equilibrium 
state, cr cx exp(—/?o-77eff), or from the observed hyper¬ 
plane E, respectively. 

Mirroring the indeterminacy of the effective Liouvil¬ 
lian, the entropy production rate 77 , too, is only deter¬ 
mined up to a multiplicative constant. When the Liou¬ 
villian is rescaled by a constant factor c, £eff 
then so is the entropy production rate, 77 —>■ cq. Up 
to this undetermined multiplicative constant, the above 
condition, Eq. (74), specifies 77 uniquely. This enables 


one to convert the proxy field W to the true vector field 
V. Erom there, the reconstruction of the dynamics can 
proceed as before. 


V. EXAMPLE 

In this section, I apply the general framework to a 
simple example. In fact, the example is so simple that 
its dynamics might just as well be analyzed without 
the whole apparatus of differential geometry introduced 
above. However, it serves to illustrate the internal con¬ 
sistency of the approach, and has the advantage of being 
solvable analytically. 

To be specific, I consider an exchangeable assembly 
of qubits, possibly interacting with each other and with 
some unknown environment. (A physical realization 
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might be a para- or ferromagnet.) The three Pauli op¬ 
erators ax, ay, and a^ constitute the relevant observ¬ 
ables; the associated manifold of macrostates is thus the 
Bloch sphere. Measurements on samples taken from the 
assembly yield the expectation values of the relevant ob¬ 
servables, X := (ax), y ■= (ay), and z := (cr^). Driven 
by some Markovian dynamics whose form is yet to be 
inferred, these expectation values change, tracing out or¬ 
bits in the Bloch sphere. The geometric shape of these 
orbits is given but —for lack of time information— their 
parametrization is not. 

For the sake of concreteness, I presume that when¬ 
ever z = 0 initially, it is z = 0 on the entire orbit; so 
orbits with this initial condition are constrained to the 
hyperplane S = {/r G 7 r( 5 )|z = 0}. Within this hyper¬ 
plane, orbits are spiralling towards the equilibrium state, 
a, which I take to be the totally mixed state, cr = //2 
(Fig. I). I assume that near equilibrium the orbits have 
the shape of a logarithmic spiral, 

(j) = (j)o+ \nr, (76) 

where r, 4> are polar coordinates in E, 

X = r cos (ft , y = r sin cj). (77) 

Given only this information, is it possible to reconstruct 
the dynamical law near equilibrium? 



FIG. 2. Two-dimensional section (z = 0) of the Bloch sphere. 
Near equilibrium, a = 7/2, orbits have the form of a logarith¬ 
mic spiral. 

Near equilibrium the entropy function is approxi¬ 
mately quadratic. 


(78) 


With entropy as a proxy for time, the orbit can be given 
a parametric representation, 

r{S) = V2\n2-2S, = (j)o+^\n{2ln2-2S), (79) 


which gives rise to the proxy field 

r dr dcf) 


(80) 


Mapping vectors to one-forms with the help of the metric 
tensor, given near equilibrium by 


C « dr (g) dr -h r^dcj) (g) d(j). 


(81) 


yields the one-form 

oj = -dr + d(j). (82) 

r 

From the orientation of the hyperplane S one readily 
concludes that oc Exploiting moreover that near 

equilibrium it is « —x (and likewise for A^), one finds 

{iCesGa)^(^p)dX^ « cr’^dcj), (83) 


where c is some undetermined multiplicative constant. 
Inserting these results into the equation for the entropy 
production rate, Eq. (74), yields 


1 

2 


dr] 


-dr + (77 
r 


cr'^)d4>. 


(84) 


which has the unique solution 


T] = cr^. 


(85) 


This entropy production rate allows one to convert the 
proxy field W to the “true” vector field V, 


V = -c 



A 

dcf) 


( 86 ) 


Einally, by Eq. (69), this yields the tensor field T gov¬ 
erning the Markovian dynamics, with the antisymmetric, 
non-dissipative part 


/_•) c \ d d d d 
r\]h® ^ 

and the symmetric, dissipative part 

tW « cC-^. 


(87) 


( 88 ) 


VI. DISCUSSION 

In the preceding sections I showed how it is possible, 
in principle, to infer a Markovian equation of motion 
— with both its non-dissipative and dissipative parts — 
from experimental data even when the latter lack infor¬ 
mation about time. To achieve this inference, I exploited 


S'«ln 2 -rV 2 . 
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the Riemannian geometry of the manifold of macrostates, 
as well as the structure of generic Markovian dynamics. 
In addition, I made a number of assumptions: 1 presumed 
that (i) energy is conserved; (ii) the experimental data 
have been taken near equilibrium (limiting thus the va¬ 
lidity of the inferred equation of motion to the vicinity 
of equilibrium, too); (iii) the dynamics is genuinely ir¬ 
reversible, S > 0; (iv) the equilibrium is canonical; and 
(v) the effective, time-independent Hamiltonian featur¬ 
ing in this canonical equilibrium state also governs the 
effective non-dissipative dynamics near equilibrium. The 
last assumption allows one to infer immediately the non- 
dissipative part of the Markovian dynamics (up to a mul¬ 
tiplicative constant); one simply identifies the effective 
Hamiltonian with the logarithm of the equilibrium state 
or, equivalently, with the observable that is conserved in 
the hyperplane S. There remains then the task of infer¬ 
ring the dissipative part of the dynamics. This requires 
more effort but is feasible, too, with the tools provided 
above. 

The reconstruction scheme fixes time intervals only up 
to a multiplicative constant, and hence time itself only 
up to an affine transformation. This freedom ensures the 
compatibility of the reconstruction scheme with special 
relativity. Two observers, moving relative to each other, 
may see the same unparametrized orbits traced out by a 
Markovian process in the manifold of macrostates. In the 
example discussed in Sec. |Vj one observer might be in 
the rest frame of the experiment, while the other observer 
moves relative to this frame in the 0 direction; both ob¬ 
servers see the same logarithmic spiral depicted in Fig. 

Their respective reconstructions of time are allowed to 
differ by an affine transformation. This provides, in par¬ 
ticular, for the possibility that they differ by a Lorentz 
transformation, as demanded by relativity. 

The time reconstructed via the above scheme is a 
macroscopic time, for it is inferred from experimental 
data on some coarse-grained level of description. If the 
same process is observed on two different levels of de¬ 
scription, and if on both levels the process is Marko¬ 
vian, then the above procedure can be applied to both 
sets of data, yielding respective macroscopic times. One 
may envision such a situation when one observes, say, a 
fluid on either the Boltzmann or the Navier-Stokes level 
of description; or when one observes a spin system on 
different coarse-grained length scales, corresponding to 
varying block spin sizes. As long as the assumptions 
spelt out above hold on all these levels of description, 
the various macroscopic times agree (up to affine trans¬ 
formations). The inferred time is thus invariant under a 
change of macroscopic scale, provided that on all scales 
the dynamics falls within the basic structure of generic 
Markovian dynamics; and provided that, moreover, on all 
scales there exists an effective Hamiltonian for its non- 
dissipative part. In this sense, the scale invariance of 
macroscopic time presupposes renormalizability [30) . 

The results in the present paper show that, in principle, 
any Markovian process can serve as a clock. Absent an 


external reference clock, one may take an arbitrary, gen¬ 
uinely irreversible Markovian process, apply the above 
reconstruction scheme, and with the help of the inferred 
equation of motion, mark points on the orbits such that 
a segment between two successive points corresponds to 
some fixed time interval. When the Markovian process is 
repeated, the system evolves along one of the orbits; and 
whenever its macrostate reaches one of the designated 
points, the clock “ticks.” In contrast to everyday clocks, 
which are based on reversible, periodic motion, such a 
clock would be based on the irreversibility of the process 
involved — similar to, e.g., radiocarbon dating. 

Conceptually, one may wonder why a reconstruction 
of time is necessary at all, and why one does not sim¬ 
ply stick to entropy as a proxy for time. In fact, any 
strictly monotonic function of time — like entropy in a 
genuinely irreversible process — serves the purpose of or¬ 
dering events; and if such an alternative parameter were 
adopted universally, it would still be possible to keep ap¬ 
pointments. What distinguishes physical time from all 
other conceivable parameters is that when expressed as a 
function of time, equations of motion become particularly 
simple. As Henri Poincare put it succinctly |31j . 

Time should be so defined that the equations 
of mechanics may be as simple as possible. In 
other words, there is not one way of measur¬ 
ing time more true than another; that which 
is generally adopted is only more convenient. 

Determining physical time is thus tantamount to finding 
the simplest possible parametrization for the widest pos¬ 
sible range of processes; where “simplest” and “widest” 
must be suitably operationalized. On a speculative 
note, the reconstruction scheme expounded here may 
be viewed as a proposal for such an operationalization: 
“Simplicity” would be embodied in the assumptions sum¬ 
marized in the opening paragraph of this discussion; 
whereas “widest” would signify the scale invariance al¬ 
luded to above. Physical time would then be singled out 
as being a fixed point under changes of scale. 

I see several avenues for further research. First, it 
will be important to apply the reconstruction scheme to 
real or simulated experimental data that pertain to more 
complex, real-world problems. Secondly, the framework 
should be extended to deal with processes further away 
from equilibrium. And finally, it might be worthwhile 
to explore in more detail the conceptual issues raised in 
this discussion, i.e., the relativistic covariance of the pro¬ 
cedure, its invariance under a change of scale, and the 
possible import of these findings on the fundamental def¬ 
inition and meaning of time. 
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